Ecology and Evolution 

High-throughput SNP-genotyping analysis of the 
relationships among Ponto-Caspian sturgeon species 

Sergey M. Rastorguev 1,2 , Artem V. Nedoluzhko 1 , Alexander M. Mazur 3 , Natalia M. Gruzdeva 1 , 
Alexander A. Volkov 2 , Anna E. Barmintseva 2 , Nikolai S. Mugue 2 & Egor B. Prokhortchouk 1,3 

'National Research Center "Kurchatov Institute", Kurchatov sq. 1, 123182 Moscow, Russia 

2 Russian Federal Research Institute of Fisheries and Oceanography (VNIRO), V. Krasnoselskaya str. 17, 107140 Moscow, Russia 
3 Center of Bioengineering, Russian Academy of Sciences, 60-letiya Oktyabrya av. 7-1, 1 17312 Moscow, Russia 



Keywords 

Genotyping, glaciation, population genetics, 
SNP, species relationships, sturgeon. 

Correspondence 

Sergey M. Rastorguev, National Research 
Center "Kurchatov Institute," Laboratory of 
Genome Analysis, NBIC Center, Kurchatov 
sq. 1, 123182 Moscow, Russia. Tel: 
+7-(499)-1 967937; Fax: +7-(499)-1 967723; 
E-mail: rastorgueff@gmail.com 

Funding Information 

No funding information provided. 



Received: 28 March 2013; Revised: 30 May 
2013; Accepted: 31 May 2013 

Ecology and Evolution 2013; 3(8): 2612- 
2618 



Abstract 

Legally certified sturgeon fisheries require population protection and conserva- 
tion methods, including DNA tests to identify the source of valuable sturgeon 
roe. However, the available genetic data are insufficient to distinguish between 
different sturgeon populations, and are even unable to distinguish between 
some species. We performed high-throughput single-nucleotide polymorphism 
(SNP)-genotyping analysis on different populations of Russian (Acipenser 
gueldenstaedtii), Persian (A. persicus), and Siberian {A. baerii) sturgeon species 
from the Caspian Sea region (Volga and Ural Rivers), the Azov Sea, and two 
Siberian rivers. We found that Russian sturgeons from the Volga and Ural Riv- 
ers were essentially indistinguishable, but they differed from Russian sturgeons 
in the Azov Sea, and from Persian and Siberian sturgeons. We identified eight 
SNPs that were sufficient to distinguish these sturgeon populations with 80% 
confidence, and allowed the development of markers to distinguish sturgeon 
species. Finally, on the basis of our SNP data, we propose that the A. baerii-Vike 
mitochondrial DNA found in some Russian sturgeons from the Caspian Sea 
arose via an introgression event during the Pleistocene glaciation. 
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Introduction 

The Russian sturgeon (Acipenser gueldenstaedtii) is one of 
the most commercially valuable fish species in the Cas- 
pian, Azov, and Black Seas due to its role as a major 
caviar producer (Vlasenko et al. 1989a,b; Birstein and 
Bemis 1997) and important aquaculture species (Fig. 1). 
Three other species are closely related to A. gue- 
ldenstaedtii: the Persian sturgeon (A. persicus), the Adriat- 
ic sturgeon (A. naccarii), and the Siberian sturgeon [A. 
baerii), the least abundant and hence most valuable spe- 
cies. Together, these four species make up the Ponto- 
Caspian clade of sturgeons. However, several difficulties 
have complicated efforts to distinguish between the vari- 
ous species in this clade. Acipenser persicus inhabits the 
Caspian Sea and cannot be readily distinguished from A. 
gueldenstaedtii by mitochondrial markers (Birstein et al. 
2000; Birstein and Doukakis 2001), or by morphological 
or anatomical differences (Vlasenko et al. 1989a,b; Ruban 



et al. 2011). Thus, the taxonomic rank of A. persicus is 
disputed, with some authors considering it to be a sepa- 
rate species (Luk'yaneko et al. 1974; Putilina 1983) and 
others regarding it to be a subspecies of A. gueldenstaedtii 
(Berg 1934; Birstein and Bemis 1997; Birstein et al. 2000; 
Birstein and Doukakis 2001; Ruban et al. 2011). Although 
the Adriatic sturgeon uniquely inhabits the northern 
Adriatic Sea, it also morphologically resembles A. gue- 
ldenstaedtii (Tortonese 1989). Unlike its anadromous rela- 
tives, the Siberian sturgeon is a freshwater species. Its 
geographical distribution includes Siberian Rivers and 
Lake Baikal, and does not overlap with Russian sturgeon. 
All these species are highly endangered, and are under the 
control and protection of the Convention for Interna- 
tional Trade in Endangered Species (CITES). 

The unusual polyphyly of Russian sturgeons' mitochon- 
drial DNA further complicates their unique identification. 
Roughly a third of Russian sturgeon specimens from the 
Caspian Sea have mitochondrial DNA that is similar to 
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Figure 1. The Russian sturgeon (Acipenser gueldenstaedtii) at the 
Astrakhan fish farm. 

that of Siberian sturgeons, which is often referred to as 
the "baerii-like mitotype" (Birstein et al. 2000; Jenneckens 
et al. 2000; Rastorguev et al. 2008). However, baerii-like 
haplotypes found in the Caspian population of A. gue- 
ldenstaedtii are different from those of A. baerii in one 
substitution in the highly variable D-loop and can be 
distinguished by PCR-based test (Mugue et al. 2008). 
Two percent of Russian sturgeons in the Azov Sea have 
the baerii-like mitotype (Timoshkina et al. 2009). In con- 
trast, no Persian sturgeons have been found to possess 
the baerii-like mitotype (Mugue et al. 2008), supporting 
that A. persicus is genetically isolated entity. Nevertheless, 
mitochondrial DNA analysis can distinguish only the 
Siberian sturgeon from among its relatives (Mugue et al. 
2008). 

There are several explanations of the "baerii-like mito- 
type" phenomenon. Some authors propose that the baerii- 
like mitotype arose from farmed A. baerii having escaped 
into the Volga River (Jenneckens et al. 2000), whereas 
others suggest that the baerii-like mitotype and the major- 
ity A. gueldenstaedtii mitotype are ancestral mitotypes for 
the Siberian, Russian, and Persian sturgeon species (Bir- 
stein et al. 2005). It is also possible that hybridization 
between Siberian and Russian sturgeons during the Pleis- 
tocene glaciation led to Siberian sturgeon mitochondrial 
DNA introgression, as considerable changes in the species 
distributions of other organisms occurred in that period 
(Hewitt 2000). 

The identification of nuclear genome markers associated 
with specific Ponto-Caspian sturgeon species has generally 
lagged behind due to the polyploidy of the sturgeon gen- 
ome. It is widely accepted that polyploidization took place 
several times within the order Acipenseriformes. First 
whole-genome duplication presumably occurred in the 
ancestral lineage and led to formation of 2n = 120-140 
("low chromosome number") group of species. Two inde- 
pendent second rounds of polyploidization were proposed: 
one in pacific clade, and another one in Atlantic clade of 
sturgeons (2n = 240-260), leading to formation of two 
phylogenetically unrelated groups of "high chromosome 



number" sturgeon species (Birstein et al. 1993). Evidence 
suggests that low chromosome number species are tetrap- 
loids with functionally diploid genome, and high chromo- 
some number species are octoploids, with two sets of 
paralogous tetrasomic loci. 

The species relationship between Russian and Persian 
sturgeons has been addressed recently by restriction site- 
associated DNA (RAD) sequencing (Baird et al. 2008), 
but the issue remains unresolved (Ogden et al. 2013). 

Meanwhile, high-throughput single-nucleotide poly- 
morphism (SNP) genotyping provides precise sample 
assignment to its origin therefore it becomes highly valu- 
able tool against illegal fishing and mislabeling. Effective- 
ness of such approach was shown for four commercial 
marine fish species (Nielsen et al. 2012). 

In this study, we performed high-throughput genotyp- 
ing analysis of individuals from several populations of 
Russian, as well as Persian, and Siberian sturgeons to find 
SNP markers for species identification. The resulting data 
should aid in sturgeon conservation efforts and help 
prevent illegal fishing. Additionally, we suggest an expla- 
nation of the baerii-like mitotype found in some Russian 
sturgeons in the Caspian Sea region. 

Material and Methods 

RNA isolation and cDNA library construction 

A several-month-old Russian sturgeon fingerling from an 
Astrakhan fish farm (Ltd Astrakhan Fish Breeding Com- 
pany, also known as "Beluga") was suspended in RNAlater 
reagent (Invitrogen/Life Technologies, Carlsbad, CA) for 
the preparation of a cDNA library. RNA extraction from 
different tissues (representing total body RNA content) 
was done with the TRIzol reagent according to the pro- 
vided RNA extraction protocol (Invitrogen). Recent inter- 
national FP6 supported project SturSNiP revealed that 
there is no fixed segregating genetic markers between these 
species, and all informative markers are frequency depen- 
dent due to sheared polymorphism (Ogden et al. 2013). 
Our approach was different from that one of SturSNiP 
because our aim was to reveal highly polymorphic loci in 
the Caspian population of Russian sturgeon; therefore, we 
collected transcriptome information from a single individ- 
ual rather than from pooled samples from different popu- 
lations of different species. There was twofold rationale 
behind taking one fingerling for transcriptome analysis; 
first, we attempted to analyze the entire mRNA repertoire, 
so we needed artificially propagated fingerling to be sacri- 
ficed and the diversity of tissues to be sampled, and sec- 
ond, the data of allelic composition from one individual 
allowed us to avoid problem of paralogous (resulted from 
previous whole-genome duplication) loci and to discrimi- 
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nate paralogous loci from true polymorphic sites (see 
below). cDNA was prepared from total RNA by a M-MLV 
reverse transcriptase approach (Schmidt and Mueller 
1999) using the MINT cDNA synthesis kit (Evrogen, cata- 
log #SK001). The resulting "SMART" cDNA was amplified 
and normalized (reducing the levels of highly abundant 
transcripts), using the duplex-specific nuclease (DSN) 
normalization method (Zhulidov et al. 2004). 

Transcriptome sequencing and mapping 

Approximately 5 /.ig of nonnormalized cDNA and 10 /ig 
of normalized cDNA were used for SOLiD library con- 
struction and sequencing using the SOLiD™ 4.0 system 
(Life Technologies) according to the standard protocol 
(SOLiD User Guide). Reads from the normalized and 
nonnormalized libraries were mapped against the 
sequence of Russian sturgeon mitochondrial DNA (NCBI 
accession number NC_012576.1), and all reads apparently 
derived from mitochondrial DNA were excluded from 
further mapping. The remaining reads were mapped 
against the sturgeon EST database from GenBank using 
the Bowtie software package (Langmead et al. 2009). An 
EST-derived reference database was created using 
sequences bearing the "Acipenser" tag in the "species 
name" field (6219 sequences in total), and was then used 
for mapping. SNP calling was performed using the SAM- 
tools package (Li et al. 2009; Li 2011). Mismatches 
between mapped reads were catalogued in VCF files, and 
heterozygous positions were then identified with a custom 
Perl script. SNP positions were sorted by coverage, and 
those that were covered by more than 100 reads were 
selected for the further analysis. 

Sturgeon samples and genotyping 

Sturgeon samples were obtained from The Russian National 
Collection of Reference Genetic Materials (RNCRGM), 
developed by the Russian CITES scientific authority on 
sturgeon fishes and maintained by VNIRO. Population 
samples of Caspian basin Russian sturgeons originated from 
the Volga and Ural Rivers (14 specimens each), with five 
additional specimens of Russian sturgeons originating from 
the Azov Sea; 28 Persian sturgeon specimens originated 
from the northern Caspian Sea; and five Siberian sturgeon 
specimens originated from the Lena and Irtysh rivers. DNA 
was isolated using the Wizard® SV 96 Genomic DNA Puri- 
fication System (Promega, Fitchburg, WI). 

Custom GoldenGate Genotyping Panel was constructed 
in Illumina Inc., San Diego, CA using loci, selected from 
SNP calling procedure described above (Table SI). The 
GoldenGate assay was performed using the Illumina 
GoldenGate platform according to the manufacturer's 



instructions and analyzed by the genotyping module of 
Illumina's Genome Studio data analysis software (Illu- 
mina Inc.). Loci for which all specimens showed the same 
genotype or genotypes could not be distinguished by this 
software, and were excluded from further analysis. 

Genotyping data analysis 

An UPGMA (unweighted pair group method with arith- 
metic mean) tree of sturgeon populations was constructed 
using the POPTREE2 software (Takezaki et al. 2010), with 
F S r distances according to Latter (1972). We used the 
Whichloci program to identify SNP loci with the highest 
statistical power for population assignments (Banks et al. 
2003). For the Whichloci analysis, we generated four simu- 
lated datasets by random resampling the SNP loci (three 
with sample size 100, one with sample size 500), and we 
then ran Whichloci on each dataset and chose loci that 
were present in the outcomes from all four simulated data- 
sets. Population structure was inferred using the STRUC- 
TURE software (Pritchard et al. 2000), with parameters 
"prior sample location information" and "allele frequency 
model correlated among populations" included in the 
admixture model. P ST statistics - the effect of subpopula- 
tions (s) compared to the total population (t) - were 
determined using the GenePop web service (http://gene- 
pop.curtin.edu.au/; Raymond and Rousset 1995; Rousset 
2008), with the "f S T and other correlation" options. 

Results and Discussion 

RNA-Seq analysis of the sturgeon 
transcriptome 

Transcriptome data for the Russian sturgeon fingerling 
were generated from approximately 250 million 50-base 
reads of a normalized cDNA library and 80 million reads of 
a nonnormalized cDNA library using the SOLiD (Life 
Technology) platform. After removing mitochondrial DNA 
reads, approximately 1.7 million reads from the normalized 
library and 3 million reads from the nonnormalized library 
mapped to the GenBank sturgeon expressed sequence tag 
(EST) database (Table 1; note that the reduced number of 
mapped reads for the normalized library was expected due 
to the removal of ribosomal and tRNA transcripts during 
the normalization process; 0.3% [821,892 of 252,186,716] 
of normalized library reads and 4.6% [3,624,935 of 
79,030,219] of control library reads were mapped on ribo- 
somal database). From 6219 totals, 4300 ESTs were covered 
by at least one read. The overall low numbers of mapped 
reads could be explained by the fact that the database con- 
tains mainly A. sinensis and A. transmontanus ESTs, which 
are distantly related to A. gueldenstaedtii. 
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SNP discovery for sturgeon genotyping and 
population identification 

Several thousand SNP candidates were identified by 
screening heterozygous positions using the following 
criteria of read sequence quality, read mapping quality, 
and quality score at the mismatch position. Additionally, 
we required that the presumptive SNPs have allele ratios 



Table 1. SOLiD read mapping and SNP calling statistics for two 
Acipenser gueldenstaedtii cDNA libraries. 





Normalized 
cDNA library 


Nonnormalized 
cDNA library 


Total reads 


252,186,716 


79,030,219 


Mapped reads 


1,787,162 


3,159,679 


Total SNPs called 


1,014,088 


755,147 


Heterozygous SNP 


50,124 


30,682 


SNPs (coverage >10) 


611,971 


415,411 


Heterozygous (coverage >10) 


45,645 


27,143 


SNPs (coverage >100) 


181,683 


138,058 


Heterozygous (coverage >1 00) 


17,765 


9812 



Read mapping was performed using the sturgeon EST database as a 
reference. SNP calling was done with default SNP caller settings (using 
SAMtools) and with minimum coverage depths of 10 and 100. 



close to 3:1 based on the following assumptions: "high 
chromosome number" sturgeons, including three species 
in our study, have two sets of paralogous tetrasomic loci 
(total eight similar genomic regions, segregated as two 
tetrasomic loci), special precaution was applied to avoid 
erroneous considering alternatively fixed paralogous loci 
as a polymorphic locus. Thus, alleles with 1:1 ratios would 
likely represent divergent paralogous loci which would not 
be polymorphic across individuals in a population. There- 
fore, we sought to identify genotypes where one tetrasomic 
set of alleles was heterozygous and the other locus (paral- 
og) was homozygous. Allelic ratio 3:1 would likely repre- 
sent situation when one paralog is fixed (i.e., AAAA) and 
the other one is polymorphic (AABB). Total 384 of the 
best-scoring mismatch positions were chosen for a Gold- 
enGate custom SNP panel (Illumina) for sturgeon geno- 
typing (Table SI). The GoldenGate assay revealed that 123 
of the 384 variable loci were successfully genotyped; the 
others appeared homozygous or heterozygous across all 
individuals (fixed at the same or alternative alleles at par- 
alogous loci) and were excluded from further analysis. 
Loci that appeared heterozygous across all individuals 
most likely reflected the presence of a paralogous locus, 
creating the appearance of two different alleles. 



Table 2. DNA sequences of the 12 SNP candidate markers used for sturgeon population assignment. 



Marker name 


DNA sequence 


EV824350.93 


TTCGATACGATAAGCCGCAATGTCTGCAGGAACACTGGCTAAACCCGCAATGCGCGGTCT[A/G] 




CTCGGTAAACGTCTACGATTCCACCTGGCTGTTGCTTTCACGCTGTCTTTGGCTGCAGCA 


ES698489.358 


ATTTrAAAGTGGARCCCTTTGTTGTTGATTCCTGGTTTAAATTTCCAGCACCTTCACAAA[A/T]TC 




CAGTAAGAAAGACATGTGAAGTAGACMSTGGCCGTGTATTTCAGATGTGCTCTGATAA 


DR977078.72 


TTGCCCAGGTCTTGAACAATATGCTATCAAGAAGTTTGCAGAGGCTTTTGAAGCCATTCC[A/G]C 




GTGCTCTTGCAGAGAACTCTGGTGTCAAGGGCAATGAACTTATCTCCAAACTGTATGCT 


DR976209.314 


GAACAGTGCACTCCACAAGAATTG I I I I I AGCCGTATAACTTCAAGAACAGAAAACTAAG[T/C] 




TGCATTATGAGTGGAGGAAGACAGATGACTGCGCTACTGTACACCCCTTTTACAACAASA 


DR977006.239 


CTCTTTCACCTTTGAGCTTCGTGATACTGGCAGATATGGATTCCTCCTTCCAGAGTCGCA[A/G]A 




TCAGGCCAACTTGCCAAGAAACAATGCTGGCAGTCAAATACATTGCCAAACATGTTCAG 


EV824380.474 


CCAAGGGCAATGTACCGCAACTTCTGGAGCCCATCCCATACGAGTTCATGGCGTAAAGCA[T/C] 




GGTCAAGCCAGAATAAAGTCCTGTTTTCAACAAAAAAAAAAAAAAAAAAAAANAAANNNA 


DR976578.395 


TTTCCCCATCAATGTGGTCGTTCAGGAAACTGGGTCTCTGGTGGAAATCAGGAATTTCTT[A/G]G 




GAGAAAAGTACATCCGCAGAGTACGCATGAGGCAGGGTGTGACCTGTAATGTTTCCCAG 


DR975536.188 


TAGAAGAGGGCGCCTGACCAAGCACACCAAGTTYGTTCGGGACATGATCCGCGAGGTGTG[T/C] 




GGCTTCGCCCCCTACGAGAGACGCGCCATGGAGCTGCTGAAGGTCTCCAAGGATAAGCGM 


DR977561.396 


GAAGTTTGCCTGCAATGGGACTGTGATTGAGCATCCAGAATATGGTGAAGTGGTTCAGCT[A/G] 




CAGGGTGATCAGCGCAAGAATATCTGCCAGTTCCTTACTGAGATTGGCTTGGCTAAGGAG 


ES698197.350 


AACAAGCCAAGGACAAGCCCCAGGAGAAACCCAAATAACATGGCAAGAGRACTAACCATC[A/G] 




CTTACAAGAGAAGTTTCAATCTCCGACACAAGCTTCCTGTCTGGGATTTCATTTTCMTTT 


DR974886.648 


TTGGTACCGTGTATCTCTCTGCTCTGGCTTTATAATGATGGGTGTCACCGTATATGAAGG[A/G]G 




GCACCTGGTAATGCASATGGTAAACCAATCAGTAATGTAAAGAATTGGATGCAATAAAG 


ES697699.316 


GAGAAGATYGACCTGAAGTTCAACCACCTCCAAGTTCGGACACGGACGCTTCCAGACCGC[T/C] 




GAAGAGAAGAAGGCGTTCATGGGACCACTCAAAAAGGACCGAATCCTCAAGGAGGAGACT 



Each polymorphic position is shown in square brackets. Sequence names consist of the GenBank accession number of the candidate EST 
sequence, with the SNP position within the EST given after the period. 
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Analysis of sturgeon population structure 

To perform a population structure analysis, genotype calls 
were filtered from the Illumina Genome Studio report 
into the STRUCTURE (Pritchard et al. 2000) and Gene- 
Pop (Raymond and Rousset 1995) programs. The Which- 
loci program (Banks et al. 2003) guided the choices of 
the most informative loci for population assignment. We 
detected few genetic differences at the 123 diagnostic 
positions between the Volga River and Ural River popula- 
tions of Russian sturgeons. Eight of the most informative 
SNP loci were sufficient to assign specimens to the Cas- 
pian Russian (Volga and Ural River) populations, Russian 
Azov Sea population, and Persian populations with 80% 
accuracy. A total of 12 loci were sufficient to make these 
population assignments with 90% accuracy (Table 2). 

Using the allele frequencies within each population for all 
123 polymorphic SNP loci (taken from the Genome Studio 
report), we constructed a population tree (Fig. 2). As 
expected, the Siberian sturgeons appeared more distant 
from the other species of the Ponto-Caspian clade. Our data 
also indicated that the Persian sturgeon was separated from 
all three samples of the Russian sturgeon, suggesting that A. 
persicus is genetically different from all A. gueldenstaedtii 
populations. These data do not support an opinion of Ru- 
ban et al. (2011) that A. persicus is an ecological morph 
within the Caspian Sea population of Russian sturgeon, but 
not a valid species. The population assignment analysis 
(STRUCTURE) also supported the genetic separation of the 
Persian sturgeon, as a homogeneous cluster consisted exclu- 



sively of A. persicus specimens (Fig. 3). The Siberian and 
Persian populations were clearly separated into two distinct 
groups, whereas the Russian sturgeons appeared as interme- 
diate and shared population structure elements of both the 
Siberian and Persian sturgeons. This is likely due to nature 
of our SNP loci when polymorphic loci were obtained from 
Russian sturgeon specimen. 



79 
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Figure 2. UPGMA tree showing the relationships of sturgeon 
populations (POPTREE2; Takezaki et al. 2010). Genetic distances (F ST ) 
between the populations (Latter 1972) were estimated using allele 
frequency data for 123 loci. Bootstrap values are indicated next to 
corresponding nodes. GUE_azov, Acipenser gueldenstaedtii from the 
Azov Sea; GUE_ur, A. gueldenstaedtii from the Ural River; GUE_vol, 
A. gueldenstaedtii from the Volga River; PERS, A. persicus; BAE, A. 
baerii. (The Ural and Volga River sturgeons are considered to be 
representative of Caspian Sea Russian sturgeons.) 
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Figure 3. Population structure analysis of 
various sturgeon samples using the 
STRUCTURE software. Three runs of 
STRUCTURE with the cluster number (K) set to 
two, three, or four are shown. Siberian 
{Acipenser baerii) and Persian (A persicus) 
sturgeons form their own distinct clusters, 
whereas all three Russian sturgeon populations 
have mixed origins. 
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Siberian sturgeon introgression into 
Caspian Russian sturgeon population during 
the Pleistocene glaciation 

Frequency of several SNPs of the Caspian (Volga and Ural 
Rivers) populations of Russian sturgeon shows an interme- 
diate value compared to those in Siberian and Persicus pop- 
ulations, indicating that Caspian population of the Russian 
sturgeon retains traces of ancient introgression of Siberian 
sturgeon into ancestral Caspian population of the Russian 
sturgeon during the Pleistocene glaciation. These data sup- 
port the hypothesis that the baerii-like mitotype is a result 
of this ancient introgression. The climatic oscillations dur- 
ing the glaciation caused great changes in species distribu- 
tion, and considerable evidence suggests that other 
organisms dispersed to new locations (Hewitt 2000). An 
investigation of 41 North American fish species concluded 
that the Pleistocene glaciation had a great impact on fish 
population structure (Bernatchez and Wilson 1998). The 
presence of distinct mitochondrial DNA sequences within a 
species has been explained as a consequence of intraspecific 
lineage reorganization, whereby the dispersed lineages of an 
ancestral set of refugees reunite in one formerly glaciated 
region (Taberlet and Bouvet 1994; Santucci et al. 1998). 

Notably, the hypothesis of the interspecific hybridiza- 
tion is required to explain abundance of the baerii-like 
mitotype and nuclear markers in the Caspian population 
of A. gueldenstaedtii. This scenario is supported by evi- 
dence that hybridization between Russian and Siberian 
sturgeons leads to the formation of fertile hybrids (Jen- 
neckens et al. 2000), which are widely used in aquacul- 
ture. Considering that the glaciation allowed the Siberian 
sturgeons into the Caspian Sea, there may have been 
nothing to prevent their hybridization with Caspian stur- 
geons. Previously published data that "baerii-\\ke" mito- 
types are present in the Caspian Russian sturgeon 
population and absent from the Persian population are 
consistent with this notion (Mugue et al. 2008). 

Conclusions 

In this study, we discovered SNP markers for distinguish- 
ing between Caspian and Azov Sea populations of Russian 
sturgeons, as well as between Russian, Persian, and 
Siberian sturgeons. These data allow the development of a 
genetic test that can be used for law enforcement in stur- 
geon fishery, aquaculture industry, and trading, as well as 
in conservation efforts for endangered sturgeon species. 
Finally, our SNP analysis supported the idea that the 
mitochondrial DNA polyphyly of Russian sturgeons in 
the Caspian Sea resulted from an introgression event 
during the Pleistocene glaciation. 
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